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ABSTRACT 

We present two supernovae (SNe) discovered with the Hubble Space Telescope (HST) in the Cosmic 
Assembly Near-infrared Deep Extragalactic Legacy Survey, an HST multi-cycle treasury program. We 
classify both objects as Type la SNe and find redshifts of z = 1.80 ± 0.02 and 2.26to;io) the latter of 
which is the highest redshift Type la SN yet seen. Using light curve fitting we determine luminosity 
distances and find that both objects are consistent with a standard ACDM cosmological model. 

These SNe were observed using the HST Wide Field Camera 3 infrared detector, with imaging in 
both wide- and medium-band filters. We demonstrate that the classification and redshift estimates are 
significantly improved by the inclusion of single-epoch medium-band observations. This medium-band 
imaging approximates a very low resolution spectrum (A/AA <100) which can isolate broad spectral 
absorption features that differentiate Type la SNe from their most common core collapse cousins. 

This medium-band method is also insensitive to dust extinction and (unlike grism spectroscopy) it 
is not affected by contamination from the SN host galaxy or other nearby sources. As such, it can 
provide a more efficient - though less precise - alternative to IR spectroscopy for high-z SNe. 

Subject headings: cosmology: observations, methods: observational, supernovae: general 


1. INTRODUCTION 

Type la supernovae (SNe) at redshifts z > 1 are a 
valuable tool for measuring the expansion history of the 
universe and for testing models of Type la SN progenitor 
systems. Measurements of the Type la SN rate at z > 1 
have been used to constrain Type la SN proge nitor mod¬ 
els through the delay time distribution test (Poznanski 


et al.|2007l|Graur et al.|20111|Dahlen et al.|2UU4[|Strofgey 

fc: Riess|2006 Dahfen et af.|2008[) and recent extensions of 
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the sample to z ^ 2 are providing unique leverage on the 
fraction of Type la SN progenitors that explode promptly 


after formation (Graur et al. 2014 Rodney et al. 2014 


hereafter |R141. 

A subset of these high-z Type la SNe with well sampled 
light curves has been used to extend luminosity distance 
measurements into an era when the universe was decel- 


erating. A notable early exa mple is SN 1997ff (Gilliland 
et al.|[l999 Riess et al. 20011, which held the title of the 


farthest kno wn 'Type la SM tor over a decade (until Jones 
et al. I 20131. By reaching into the early universe, SiN 
1997n and other high-z Type la SNe provided a valuable 
check against systematic effect s that could mimi c dark 
energy signatures at lower z dRiess e t al.||2000l |2001 
Tonry et al. 2003 Riess et al. 20046| l As the sample 


of high-z Type la SiNe has grown, it has also delivered 
limits on the evolution o f the equation of state of the uni- 
verse over cosmic time dRiess et al. 2007 Suzuki et al. 


20121 . 

Discovering Type la SNe and measuring their rest- 
frame optical light curves at z > 1 requires deep infrared 
imaging. This imaging data alone has been widely used 
for measuring Type la SN rates to test SN progenitor 
models, aided by an increasingly r igorous suite of photo¬ 
metric classification software (e.g.,|Poznanski et al.|2007[ 
' Rodney &: Tonry||2009| [Kessler et al.||2010[ |Sako et al'| 


2011D . 


However, SN cosmology is still almost exclusively lim- 
ited to samples with spectroscopic cla ssificatio ns (e.g.. 
Rest et ail 120141 [Betoule et al.||20T4| but see jCamp-j 


beil et al.||2013|). Spectroscopic observations of SNe in 
the range 1 < z < 1.5 with ground-b ased telescopes 


2005 


ageat 


Morokuma et al. 

2010 

). The task is more man- 

ie from space, and 

slitless grism spectroscopy with 
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the Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) has been used to su ccessfully classify 16 
Type la SNe in this redshift range (Riess et al.||2004a|b 
Riess et al. |2007t Graur et al. 2014'R14). However, as 
SiNe discoveries are pushed to redshitts z > 1.5, such di¬ 
rect spectroscopic SN cla ssifications become increasingly 
costly and less definitive (Rodney et al. 2012 Rubin et al. 


2013 


Jones et ^[2013 ). 


In cases where the SN spectrum cannot independently 
define the SN class and redshift, we may make use of pho¬ 
tometric classifications and host galaxy redshifts. How¬ 
ever, biases in the derived redshift and distance are more 
likely in such cases, due to errors in host galaxy as¬ 
sociation or when a spectroscopic redshift for the host 
galaxy is unavailable. Host galaxy redshifts are espe¬ 
cially problematic for SNe in the so-called redshift desert 
(1.4 < z < 2.5), which is precisely the current frontier 
for high-z SN searches with HST. Newer infrared spec¬ 
trometers on 8 -m class telescopes are able to define pre¬ 
cise spectroscopic redshifts for some SN host galaxies in 



Amanullah et al. 

2011 

Frederiksen et al. 

2012 

2014| 

. These biases can also be reduced - though not 


clusters, where SNe appearing in early-type cluster mem- 
ber galaxies have less ambiguity in both class and redshift 
( Dawson et alj2009 ). 

in Section we present an alternative approach for 
determining the classification and redshift of SNe at 
z > 1.4. This method builds on the “SN Cross Correla¬ 


tion” imaging technique of Scolnic et al. (2009 hereafter 
S09), extending it to higher redshitts by using medium- 
band imaging from the HST Wide Field Camera 3 in¬ 
frared (WFC3-1R) detector. Unlike photometric con¬ 
straints from pure broad-band imaging, this approach 
is insensitive to reddening, and is able to test for spe¬ 
cific individual absorption features in the SN spectral 
energy distribution (SED). This medium-band method 
is more efficient than slitless spectroscopy with an HST 
grism, and is also largely unaffected by contamination 
from nearW bright sources. 

Section describes the discovery and observation of 
two Type la SNe at z ~ 2, found in the SN search com- 


tragalactic Legacy Survey (CANDELS 

Koekemoer et al. 

2011 
the T 

Grogin et al.||2011 
St medium-band i 

). In Sections 
making: techni 

Que ton 

we apply 
elD define 


the classification and redshift of each SN. A third high- 
redshift SN, classified as a Core Collapse SN (CCSN), is 
presented in the Appendix. The value of the medium- 
band imaging method is examined in Section and in 
Section we put these SN discoveries into the context of 
current Type la SN cosmology. Finally, our results and 
conclusions are briefly summarized in Section 

2. PSEUDO-SPECTROSCOPY WITH 
MEDIUM BAND FILTERS 

The SED of a Type la SN near maximum light is 
dominated by a series of broad absorption features inter¬ 
spersed with pseudo-continuum peaks, creating a roughly 
sinusoidal pattern over the rest-frame wa velength range 

introduced the 


3500-6500 A (e.g., |Filippe^|l997 |. 


S09 


specialized Superb]ovA Cross Correlation (SNACC) fil¬ 
ters, a set of custom optical filters designed to take ad- 
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Figure 1. The three primary medium band filters for the HST 
WFC3-IR camera, and their complementary wide band filters. The 
transmission function of each medium-band filter is shown as a 
filled curve, and the accompanying wide-band filter as a solid black 
line. These three pairs define the medium-minus-wide pseudo¬ 
colors used throughout this work. 


vantage of this crudely periodic SED shape. Each filter 
has 3 pass-bands (“teeth”) that approximately matc h the 
width and spacing of the SNIa spectral features. |S09 


showed that SNACC filter observations combined with 
traditional broad band imaging from the Pan-STARRSl 
SN survey could deliver a Type la SN sample with 98% 
purity and redshift precision of = 0 . 01 , all with an 
efficiency up to ^30 times better than traditional spec¬ 
troscopic follow-up. In addition to the efficiency gains, 
this comb filter approach is largely insensitive to dust ex¬ 
tinction, and can easily employ image subtraction tech¬ 
niques to minimize contamination of the SN signal by 
the underlying host galaxy light. _ 

The SNACC filters employed by |S09| were optimized 
for SNe in the redshift range 0.2 < z < 0.8, in or¬ 
der to target SNe discovered in ground-based surveys. 
To extend this technique to higher redshift SNe, we 
must turn to infrared (IR) wavelengths. The HST 
WFC3-IR camera includes a set of medium-band fil¬ 
ters (F098M,F127M,F139M,F153M) which can collec¬ 
tively approximate the teeth of the SNACC comb filters. 
These four medium-band filters on the WFC3-IR detec¬ 
tor were designed as grism reference bands or to pick out 
molecular absorption features in extrasolar planets (e.g., 
H 2 O, CH 4 , NH 3 ). The F098M filter spans 157nm, while 
the three redder bands have a width of ~ 65 nm. By 
comparison, the wide-band IR filters such as F125W (J) 
and F160W (77) cover ^ 275 nm each. As shown in 
Figure each of the medium-band filters has a corre¬ 
sponding wide-band filter that completely contains the 
medium-band filter wavelengths. Throughout this work 
we use medium-minus-wide pseudo-colors, subtracting 
the magnitude of a source in the wide band from the 
medium-band magnitude (e.g., F12 7M-F 125W). This is 
consistent with the convention set in lS09l for the SNACC 
filters. 

The spectral absorption features of Type la SNe span 
roughly 30 nm in the rest frame, so for a high-redshift SN 
at 1.5 < z < 2.5 the observer finds them broadened to 
widths of 75 — 105 nm. As shown on the left side of Fig¬ 
ure this means that the medium-band filter widths are 
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Figure 2. Demonstration of the method for determining a Typ e la SN redshift from medium—wide band pseudo-colors. Three panels at 
left show the SED of a template Type la SN (|Hsiao et al.|2007j at peak brightness with no extinction. From top to bottom these panels 
show the SED at redshifts of 1.8, 2.0, and 2.2. Uveriaid on eacn SED are the transmission curves for the WFC3-IR medium band filters. 
As the redshift increases, broad absorption features are moved in and out of the band passes, causing cyclic changes in the observed flux 
through each band, relative to the encompassing broad bands (not shown). In two panels at right we plot the medium—wide pseudo-colors, 
mapping the results of changing the redshift from z = 1.8 to 2.2. 


well suited for isolating the spectral features character¬ 
istic of Type la SNe in this redshift range. Observations 
with traditional wide-band filters must accompany each 
medium-band observation in order to provide a contin¬ 
uum reference point. The medium-band pseudo-colors 
(e.g. F139M-F140W) indicate whether the medium band 
filter intersects an absorption feature (resulting in a more 
positive pseudo-color) or a brighter continuum region 
(more negative pseudo-color). As redshift increases, the 
absorption features migrate in and out of the medium 
band regions, so pseudo-colors move from positive to neg¬ 
ative and back again. In a plot comparing two pseudo¬ 
colors, as shown in Figure a Type la SN will trace 
out a circle as it is moved through redshift space. For 
redshifts 1.5 < z < 2.5, measuring three of these pseudo¬ 
colors near peak brightness to better than 0.1 mag can 
yield a redshift estimate with precision approaching 2% 
in Az/(1 + z), comparable to the best ava ilable photo¬ 


metric redshifts from galaxy SED fitting (Hildebrandt 
et al.|[2Mo{ Dahlen et al.||2013 |. 


The medium-band pseudo-colors are also useful for 
classification of SNe at z > 1.5. The spectra of core col¬ 
lapse SNe (CC SNe) do not have the same semi-periodic 
absorption features as in Type la SNe. CC SN SEDs are 
also much more heterogeneous, so the space that they 
populate in a pseudo-color versus pseudo-color plot is 
not a circle, but an amorphous cloud. As shown in Fig¬ 
ure!^ the most common sub-class of CC SNe (Type II) is 
centered near the origin of the pseudo-color space, so the 
location of a SN in this space can be a useful indicator 
of its class. 

The use of medium-band pseudo-colors has several 


qualitative features that make it attractive for SN clas¬ 
sification. First, by using relative flux instead of ab¬ 
solute flux measures (i.e. colors instead of magnitudes), 
this method does not require any assumptions about the 
underlying cosmology or luminosity function of the high- 
redshift SN populations. Second, the wide- and medium- 
band filters in each color pair occupy the same portion 
of the SED, so this method is insensitive to the over¬ 
all color of the SED. Therefore, as shown in Eigure 
the classification is largely unaffected by SN host galaxy 
extinction: medium—wide pseudo-colors are changed by 
less than 0.01 mag per 1 mag of extinction. 

Medium-band imaging is also sub stantially more efS - 
cient than spectroscopy. As shown in Jones et al. (|2013 ), 
observing a typical Type la SN at z > 1.5 with the 
WEC3-IR G141 grism requires at least 10 HST orbits 
(25 ks) to get sufficient signal-to-noise (S/N) for discrim¬ 
inating a Type la SN spectrum from that of a Type II 
SN. To rule out possible mis-classification of Type Ib or 
Ic SNe as Type la SN requires considerably more: as 
many as ~70 orbits at z ~ 2. HST grism observations 
are also often limited by systematic uncertainties due to 
contamination from the SN host galaxy and other nearby 
sources that overlap the SN spectrum. By contrast, in 
this medium band imaging method we use spectral ele¬ 
ments that have ~15% better total throughput at any 
given wavelength than the G141 grism. 

With medium-band imaging we can also utilize exist¬ 
ing wide-band template images (already needed for mea¬ 
suring the wide-band light curve) to subtract out con¬ 
taminating background light very cleanly. That is, we 
get the advantages of host galaxy subtraction without 
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Figure 3. Demonstration of the segregation of SN populations in medium—wide pseudo-color space. Contours show the colors measured 
from 5000 simulated SN of each class, evenly sampling the redshift range z=[ 1.8.2. 2l. with Monte Carlo sampling of light curve templates 
and parameters - including luminosity and host galaxy extinction - following |R14| Solid line contours enclose 68% of the population and 
filled contours enclose 95% for each class, with red for Type la SNe, green for 't ype II, and yellow for Type Ib/c. The separation of the 
Type la contours from the CC SN sub-types means that these pseudo-color plots can be useful aids for classifying observed high- 2 : SNe. 
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wavelength (/im) 

Figure 4. The medium-band pseudo-colors are largely insensitive to extinction. Black curves show the SED of a Type la SN template 
at z=2.0 and peak brightness j Hsiao et al.||2007ll. Three curves show the effect of host galaxy extinctions from Ay=0.0 to 2.0, using 
an extinction law with Ry=3.1 l|uarclelli et al!]|i989| |0’Donnell||1994|l. The SEDs have been normalized at 1 fim. Although the overall 
SED is substantially reddened, tne pseudo-colors change by only '^'ZVo, as given in the legend in AB mag colors for F127M-F125W and 
F153M-F160W. 
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requiring an additional post-SN epoch of medium-band 
imaging. To do so, we make the justifiable assumption 
that the host galaxy does not have any broad spectral 
features that match the width of the medium band fil¬ 
ters. With that assumption, we can create a SN-free 
medium-band template image by simply scaling the flux 
of the wide-band template using a ratio of the medium- 
and wide-band filter transmissions. 

In the following sections we will examine two examples 
from the CANDELS SN program where medium-band 
imaging is employed at a total cost of only 5-10 HST 
orbits (~10-20 ks) across all medium-band filters. 


3. SN DISCOVERY AND OBSERVATIONS 

We now turn from the theory of this technique to the 
actual practice, as applied to a pair of SNe discovered 
in the CANDELS multi-cycle treasury program on HST. 
The SN component of the CANDELS program was op¬ 
erated jointly alongside the Cluster L easing And Su¬ 


perno va Survey with Hubble (CLASH, ||Postman et al. 
2012 |, and the two programs together discovered and 


characterized ~100 SNe over 3 years ([Rodney et al.|2012 
Jones et al.| |2013| |R14[ |Graur et al.||2l)14| [Hatel et al. 
w-up observ; ' ~ 


2014p Follow-up observations with hlST were executed 
through the CANDELS/CLASH SN program (PLRiess, 
HST Program ID:12099). The SNe in this work were 
first presented in |R14[ where we measured the volumet¬ 
ric Type la SN rate as a function of redshift from the full 
CANDELS SN sample. 

SN GNDlSSto (Section and GND12Gol (Sec¬ 
tion!^ were discovered in IR imaging on the GANDELS 
GOODS-S and GOODS-N fields, respectively. In both 
cases these SN were immediately identified as likely high- 
redshift sources based on preliminary host galaxy photo- 
z measurements and their observed brightnesses. These 
SNe were also initially classified as Type la candidates, 
based on their rest-frame UV-optical color deficit (as in 


Riess et al.|2004a Rodney et al.|2012 Jones et al.|2013 ). 
No dehnitive spectroscopic redshitts were available for ei- 
ther of these objects at the time of discovery, and they 
were both quite faint, so spectroscopic follow-up with the 
HST grisms was deemed too risky. Instead, we applied 
the medium band imaging technique, described in detail 
below. 

Table [l] summarizes the discovery dates, coordin ates, 
and host galaxy redshifts of all three objects. In |R14| 
we provided classifications and redshifts based only on 
the broad-band SN light curves and host information. 
Here we will refine both class and redshift with the addi¬ 
tion of medium-band IR data and improved photometry. 
We also extend the analysis to determine luminosity dis¬ 
tances for the two objects classified as Type la, in order 
to compare them to the larger SN samples at z < 1. 


3.1. Data Processing and Photometry 

For each SN, we collected from the HST archive all 
available imaging from the Wide Field Gamera 3 (WFG3- 
IR and WFC3-UVIS) and the Advanced Gamera for Sur¬ 
veys Wide Field Gamera (ACS-WFG). These data were 
then processed using a custom pipeline sndrizpip^^ 
that sorts the HST images into epochs, then combines 

sndrizpipe vl.2 https://github.com/srodney/sndrizpipe 


and registers them using the DrizzlePac soft ware from 
the Space T elescope Science Institute (STScI) (|Fruchter| 
et al. 2010). For each filter, a template image was con- 
structed from all epochs that were at least 30 days before 
the SN appeared or 90 days after the last 3tT detection. 
The template images were directly subtracted from ev¬ 
ery epoch containing SN light to get difference images. 
It was not necessary to convolve with the point spread 


funct ion (PSF) for optimal subtraction (Alard & Lupton 
1998), owing to the excellent stability of the HS'i’ HSh'. 


"Fluxes and AB magnitudes for GND13Sto and 
GND12Gol are given in Tables and respectively. In 
the WFC3-IR difference images we measured the SN flux 
using PSF-fitting photometry, with a PSF model con¬ 
structed from the G2V standard star P330E. Flux un¬ 
certainties were estimated by separately planting 10,000 
synthetic SNe into the difference image and fitting a 
Gaussian to the distribution of recovered fluxes. 

For both of these SNe we have no detections above 
3cr in the ACS-WFG or WFC3-UVIS imaging. In those 
bluer bands (corresponding to the rest-frame UV for 
these SNe) we used the hstphot package^ to collect 
aperture photometry in a 0.4" aperture, witnzero points 
and aperture corrections as published by STScIp^ 

3.2. Galaxy Photometry, Host Association 
and SED Fitting 

To evaluate galaxies in the vicinity of our SNe as can¬ 
didate hosts and potential gravitational lenses, we use a 
photometry catalog including data from all available op¬ 
tical and I R bands. All photometry was collected using 
SExtractor (iBertin fc Arnouts |1996 ) for source detection 
and TFIT (Laidler et al. |2007|) tor m a tched -aper ture pho¬ 
tomet r y, as detailed injCuo et al. (2013) and Fontana 
et al.l (2014). For CAN DELS data where the SN was 
detected (i.e. the WFC3 IR bands) we measured pho¬ 
tometry from the SN-free template images instead of the 
full-depth CANDELS mosaics. 

To determine which of the nearby galaxies (if any) 
could plausibly be associated with e ach SN, we follow 
the procedure of|Sullivan et al.|(|2006|), which is designed 
to infer the most likely SN host galaxy based on angu¬ 
lar separation and galaxy s ize. Using SE xtractor on the 
F160W images, we define a Kron (1980) ellipse centered 


at position Xgai, ygai, with semimajor axis r^, semiminor 
axis rB and position angle 9. The separation from the 
galaxy to any given position {x, y) is quantified by the 
elliptical radius, R, defined as 


R — Cxx^r A ^yyVr A 


( 1 ) 


where Xr = x - Xgai, Ur = y - Vgai, Cxx = cos‘^(9)/r^ + 
sin^(0)/r^, Cyy = sin^(d)/r^ -|- cos^(d)/r^, and Cxy = 
2cos(d) sin(0)(l/r^ — l/r^). An ellipse with radius 
R = 2.5 will typically contain ~90% o f the object’s light 
( Infante|1987[ Graham & Driver 2005[). For c onsideration 


as a possible SN host galaxy, Sullivan et al. 


the requirement that the SN position must 


(2006) used 
ie within an 


https://github.coin/srodney/hstphot ba sed on a transla- 
tion of the IDL AstroLib photometry routines l|Landsman||1993|l 
into Python. 

For WFC3: http://www.stsci.edu/hst/wfc3/phot_zp_lbn 
for ACS: http://www.stsci.edu/hst/acs/analysis/zeropoints 
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Table 1 

SN Discovery Data 


Supernova 

SN Coordinates^ 

Host Galaxy Coordinates^ 

Discovery Date 

Host Redshift 

Host z Source^ 

GNDlSSto 

GNDI 2 C 0 I 

GSDllBus^ 

12:37:16.77 -1-62:16:41.4 
12:36:37.58 -1-62:18:33.0 
03:32:42.78 -27:48:07.1 

12:37:16.59 -1-62:16:43.4 
12:36:37.51 -1-62:18:32.6 
03:32:42.78 -27:48:07.1 

UT 2013 Junl4.5 
UT 2012 May27.3 
UT 2011 Aug05.8 

1.80 ±0.02 

2.26 ±0.09 

-*-•-*-^-0.38 

Spitzer IRS 
GANDELS phot-z 
GANDELS phot-z 


R.A. and Decl. in J2000 ICRS coordinates. 

^ The center of the SN host galaxy in J2000 ICRS coordinates. For GNDlSSto we give coordinates of the most likely host: 
stoneA. 

^ Primary source for host galaxy redshift, independent of the SN photometry. See text for details. 

^ SN GSDllBus is discussed in the appendix. 


Table 2 

GNDlSSto Photometry 


Obs. Date 
(MJD) 

Filter 

Exp. Time 
(sec) 

Flux 

(counts/sec) 

Flux Err 
(counts/sec) 

AB Mag^ 

Mag Err 

AB Zero Point 

AZPb 

(Vega-AB) 

56403.4 

F350LP 

434 

0.0236 

0.1375 

>27.9 


26.949 

-0.156 

56457.5 

F350LP 

434 

0.1084 

0.1430 

>27.9 


26.949 

-0.156 

56490.6 

F350LP 

544 

0.1271 

0.1292 

>28.0 


26.949 

-0.156 

56513.4 

F350LP 

434 

-0.0748 

0.1548 

>27.8 


26.949 

-0.156 

56406.7 

F814W 

6708 

0.0864 

0.0675 

28.606 

0.849 

25.947 

-0.424 

56458.3 

F814W 

2236 

0.1512 

0.0876 

27.998 

0.629 

25.947 

-0.424 

56410.7 

F850LP 

1946 

-0.1063 

0.0917 

>26.3 


24.857 

-0.519 

56403.4 

F125W 

1156 

0.2437 

0.1155 

27.763 

0.515 

26.230 

-0.901 

56457.5 

F125W 

1156 

1.1137 

0.1564 

26.113 

0.152 

26.230 

-0.901 

56467.9 

F125W 

1256 

2.1157 

0.1614 

25.416 

0.083 

26.230 

-0.901 

56490.6 

F125W 

1006 

2.3250 

0.1907 

25.314 

0.089 

26.230 

-0.901 

56513.4 

F125W 

1156 

1.2606 

0.1651 

25.979 

0.142 

26.230 

-0.901 

56531.1 

F125W 

503 

0.4965 

0.1907 

26.990 

0.417 

26.230 

-0.901 

56474.4 

F139M 

5023 

0.4202 

0.0791 

25.420 

0.204 

24.479 

-1.076 

56474.5 

F140W 

1356 

2.7573 

0.3110 

25.351 

0.122 

26.452 

-1.076 

56474.9 

F153M 

5023 

0.4106 

0.0621 

25.429 

0.164 

24.463 

-1.254 

56403.4 

F160W 

1206 

0.1636 

0.1493 

27.911 

0.991 

25.946 

-1.251 

56457.5 

F160W 

1206 

1.0590 

0.1611 

25.884 

0.165 

25.946 

-1.251 

56467.9 

F160W 

1206 

1.2366 

0.1543 

25.715 

0.136 

25.946 

-1.251 

56474.5 

F160W 

1206 

1.9397 

0.1501 

25.227 

0.084 

25.946 

-1.251 

56490.6 

F160W 

1006 

1.5039 

0.1728 

25.503 

0.125 

25.946 

-1.251 

56513.4 

F160W 

1206 

1.3238 

0.1740 

25.641 

0.143 

25.946 

-1.251 

56531.1 

F160W 

1609 

0.9151 

0.1628 

26.042 

0.193 

25.946 

-1.251 

56551.2 

F160W 

2462 

0.7125 

0.1277 

26.314 

0.195 

25.946 

-1.251 


^ For flux values of less than Icr significance we report the magnitude as a Scr upper limit 
^ Zero point difference: the magnitude shift for conversion from AB to Vega magnitude units. 


ellipse with i? = 5, and found that ~93% of SNe in the 
Supernova Legacy Survey (SNLS) had at least one de¬ 
tected host candidate meeting that criterion. Ellipse pa¬ 
rameters for the nearest galaxies to each SN are reported 
in Table Note that this approach does not make use 
of any redshift information about the SN or the nearby 
galaxies. In the following sections we will also bring to 
bear the available spectroscopic and photometric redshift 
evidence. 

To determine the photo- 2 ; and estimate physical pa¬ 
rameters, we fit each galaxy’s SED using the LePhare 


code for SEP matching via minimization (Arnouts 


et al. 

1999 

plate 

ibrarv 


(Pickles 1998), the LePhare q uasar templates, and th e 
BCU3 models for galaxy SEDs (Bruzual fc Char lot|2003 ). 
A Chabrier initial mass function is assumed ([Chabrier 


2003). Eurther details o f the SED-fitting process are 


given in Hemmati et al. (2014). Note, however, that 


Hemmati et al. (2014) each galaxy is subdivided to 


kpc-scale resolution. Eor this work we use the integrated 
light of the entire galaxy. Table reports the results 
from this SED fitting, giving the ^otometric redshift, 
best-fit galaxy type, and an estimate of the total stellar 
mass. We find the photometric redshifts are in very good 
agreement with the few available spectroscopic redshifts. 
In those cases where a spectroscopic redshift exists, the 
redshift is fixed at that value for determination of the 
best-fit galaxy type and stellar mass. 


4. GND13Sto : SN STONE 
4.1. GNDlSSto Host Galaxy 

The first of our two high- 2 : Type la SN candidates is 
GND13Sto, □ discovered in exposures taken on 14 June 

Nicknamed “SN Stone,” after Thomas Stone, a signatory of 
the United States Declaration of Independence as a delegate from 
Maryland 
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Table 3 

GNDI 2 C 0 I Photometry 


Obs. Date 
(MJD) 

Filter 

Exp. Time 
(sec) 

Flux 

(counts/sec) 

Flux Err 
(counts/sec) 

AB Mag“ 

Mag Err 

AB Zero Point 

AZP*’ 

(Vega-AB) 

56018.7 

F350LP 

434 

0.2079 

0.1657 

28.654 

0.865 

26.949 

-0.156 

56074.3 

F350LP 

868 

-0.2289 

0.1649 

>27.7 


26.949 

-0.156 

56129.4 

F350LP 

434 

-0.3159 

0.2412 

>27.3 


26.949 

-0.156 

56145.5 

F350LP 

400 

-0.0398 

0.1953 

>27.5 


26.949 

-0.156 

56070.0 

F606W 

982 

-0.0505 

0.1123 

>27.7 


26.493 

-0.086 

56016.3 

F814W 

2536 

0.1885 

0.1006 

27.759 

0.580 

25.947 

-0.424 

56074.6 

F814W 

3535 

0.1318 

0.0952 

28.147 

0.784 

25.947 

-0.424 

56116.5 

F814W 

7072 

-0.0889 

0.0746 

>27.6 


25.947 

-0.424 

56128.4 

F814W 

6708 

-0.0630 

0.0834 

>27.5 


25.947 

-0.424 

56180.9 

F814W 

2236 

-0.0433 

0.1186 

>27.1 


25.947 

-0.424 

56200.2 

F814W 

2396 

0.0735 

0.0479 

28.781 

0.708 

25.947 

-0.424 

56152.1 

F850LP 

1442 

-0.0647 

0.0746 

>26.5 


24.857 

-0.519 

56018.6 

F125W 

1156 

-0.0049 

0.0703 

>27.9 


26.230 

-0.901 

56074.3 

F125W 

1156 

1.2037 

0.1440 

26.029 

0.130 

26.230 

-0.901 

56084.7 

F125W 

1206 

1.1206 

0.1487 

26.106 

0.144 

26.230 

-0.901 

56103.2 

F125W 

1206 

0.8749 

0.1514 

26.375 

0.188 

26.230 

-0.901 

56129.5 

F125W 

2309 

0.2662 

0.0990 

27.667 

0.404 

26.230 

-0.901 

56145.6 

F125W 

2562 

0.2219 

0.1033 

27.864 

0.505 

26.230 

-0.901 

56183.1 

F125W 

1156 

-0.0035 

0.0708 

>27.9 


26.230 

-0.901 

56242.0 

F125W 

1156 

0.0000 

0.0763 

>27.8 


26.230 

-0.901 

56297.7 

F125W 

1156 

0.0140 

0.0742 

>27.9 


26.230 

-0.901 

56084.8 

F127M 

7485 

0.3853 

0.0447 

25.677 

0.126 

24.641 

-0.961 

56082.8 

F139M 

8291 

0.3472 

0.0527 

25.628 

0.165 

24.479 

-1.079 

56083.1 

F140W 

1706 

1.9370 

0.1847 

25.734 

0.103 

26.452 

-1.076 

56091.5 

F140W 

506 

1.8524 

0.2847 

25.783 

0.167 

26.452 

-1.076 

56104.2 

F140W 

506 

1.3515 

0.2988 

26.125 

0.240 

26.452 

-1.076 

56085.7 

F153M 

7485 

0.2942 

0.0543 

25.791 

0.201 

24.463 

-1.254 

56018.6 

F160W 

1206 

-0.0116 

0.1823 

>26.6 


25.946 

-1.251 

56074.3 

F160W 

1206 

1.3706 

0.1727 

25.604 

0.137 

25.946 

-1.251 

56084.7 

F160W 

1206 

1.0138 

0.1547 

25.931 

0.166 

25.946 

-1.251 

56103.2 

F160W 

1256 

0.9008 

0.1550 

26.059 

0.187 

25.946 

-1.251 

56129.4 

F160W 

2762 

0.4049 

0.1231 

26.928 

0.330 

25.946 

-1.251 

56145.5 

F160W 

2009 

0.5715 

0.1500 

26.553 

0.285 

25.946 

-1.251 

56183.1 

F160W 

1206 

0.2813 

0.1489 

27.323 

0.575 

25.946 

-1.251 

56242.0 

F160W 

1206 

0.0664 

0.2891 

>26.1 


25.946 

-1.251 

56297.7 

F160W 

1206 

-0.0234 

0.2204 

>26.4 


25.946 

-1.251 


^ For flux values of less than Icr significance we report the magnitude as a Scr upper limit 
^ Zero point difference: the magnitude shift for conversion from AB to Vega magnitude units. 


Table 4 

Kron Ellipse Parameters and Properties of Host Galaxy and Lens Candidates 


ID r-A (") r-B (") 0 (deg) itsN spec-z photo-z*^ logioiM/M©) 


colfaxA 

0.34 

0.25 

0.0 

1.6 


9 9c +0.09 
-0.09 

Q 98 +0.10 
y-yo _o.o6 

colfaxB 

0.64 

0.46 

89.4 

4.4 

0.9466 ± 0.0002 

'J-yo -0.04 

9 87 +9.17 
-0.14 

stoneA 

0.60 

0.26 

52.1 

3.9 

1.80 ± 0.02 

1 70 

-0.13 

11 37 +9-31 

11.0/ _o.i2 

stoneB 

0.59 

0.41 

76.8 

5.0 


cc +0.66 

-*-•00 _o.22 

10 13 +9-29 
10.10 _o.25 

stoneC 

0.66 

0.31 

51.1 

5.5 

0.5572 ± 0.0002 

t*-D0 _o.07 

Q 90 + 0.20 

-0.24 

bushA 

0.15 

0.13 

1.1 

0.0 


115 +0-19 
-0.38 

7 60 +9-29 
/ .OU _o.ig 


spec-z Ref. 


Cowie et al. 

2004 


Murphy et a 

1. 

201 

)9 


Wirth et ah 


2004 


^ Photometric redshift and 68 % confidence region. 

^ logjQ of the total stellar mass of the galaxy and 68 % confidence region, with redshift fixed to the spec-z when 
available. 
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SN GNDlSSto F160W 


12:37:16.778, +62:16:41.43 Template 



Rsn=12 

7—1 1 

^“-‘-■^-0.5 1 arcsec 



Redshift 


Figure 5. The host environment of SN GNDlSSto. The template 
F160W image shown is a combination of all exposures from before 
MJD=56400 or after MJD=56600. Ellipses are drawn at i? = 
2.5 for the five nearest host galaxy candidates, using parameters 
measured with SExtractor in the F160W band (see Equation 
Each host galaxy candidate is labeled with the R value of an ellipse 
that passes through the position of the SN, and the best available 
redshift constraint. As detailed in the text, galaxy A is the most 
probable host, while galaxies B, D and E all have broad photo- 2 : 
distributions consistent with the redshift of galaxy A. Some or all 
of these four galaxies may be members of a small group hosting SN 
GNDlSSto. 


Figure 6. Redshift constraints on GNDlSSto, assuming it is a 
Type la SN. The independent constraint from each of S different 
sources and the combined constraint are shown as probability dis¬ 
tribution functions over redshift, all normalized to a peak value of 
1.0. The composite redshift prior from host galaxy candidates A 
and B is given as a filled orange curve. The filled blue curve shows 
the marginal posterior probability from light curve fitting with only 
the wide-band data and no host galaxy redshift prior. The filled 
red curve shows the constraint from the single epoch of medium 
band data, along with the corresponding wide-band observations 
in the same epoch, but no host galaxy redshift prior. Finally, the 
solid black curve shows the composite redshift constraint derived 
from light curve fitting with all of the above information. 


2013 in the GOODS-N field. There are 22 galaxies de¬ 
tected in the CANDELS F160W imaging that lie within 
10" of the SN position, but none of these are coincident 
with the SN. 

Using the R parameter (Equation as an indication 
of the likelihood of physical association, we reduce the 
list of possible hosts to the 3 galaxies highlighted in Fig¬ 
ure which contain the SN position within an ellipse 
having R < 10. Of these, candidate A is the most 
likely, with R = 3.9, though host candidates B and C 
are close behind, wi th R = 5.0 and 5.5, respectively. 


Murphy et al. (2009) measured the redshift of galaxy A 
as z =1.82 ± 0.02, based on the IR spectrum from the 
InfraRed Spectrograph (IRS) on the Spitzer space tele- 
sc ope, and it was subsequent ly updated to z = 1.80±0.02 
Kirkpatrick et al. (|2012 ). This IRS spectrum was 


also revisited by Ilernan~Caballero| (2012), who used a 


maximum-likelihood spectral template fitting approach 
that confirmed th e solution around z = 1.8. However, 
[Hernan Caballero| erroneously associated the IRS spec- 
trum with SDSS J123716.59+621643.9 (our galaxy C) 
and thus incorrectly identified it as an outlier. 

Host candidate B is separated from the SN by i?sN = 
5.0, and has a photometric redshift of z = 1.55^g'22. 
Though quite broad, this photo-z overlaps the spectro¬ 
scopic redshift of galaxy A, z = 1.8. There are two other 
galaxies within 10" of the SN (D and E in Figure]^ that 
also have very broad photometric redshifts consistent 
with galaxy A: z = l-OlgJ for galaxy D and z = 1.7±0.6 


for E. These two low surface brightness galaxies are very 
unlikely to be directly associated with GND13Sto, as 
they are separated from the SN by i?sN > 10. Never¬ 
theless, the similarity of these photo-z estimates suggests 
that galaxies A and B may be two dominant members of 
a small galaxy group containing D and E. Given the weak 
redshift constraints, this is very speculative, but if true 
then it is possible that SN GND13Sto is either bound to 
an undetected member of that group, or belongs to the 
intracluster light. 

The only remaining galaxy near enough to be related 
to the SN is host candidate C. This is the galaxy SDSS 
J123716.59+621643.9, which has a spectros copic redshift 
of z = 0.5572 from an optical spectrum (Wirth et al. 


2004). We identify galaxy G as a foreground object be- 


cause the SN photometry is inconsistent with this red¬ 
shift: no normal SN light curve template at z = 0.56 can 
match the observed photometry with a P^r degree of 
freedom (xt) better than 3.3. 

In summary, we find that GND13Sto is most likely 
associated with either galaxy A or galaxy B. If the for¬ 
mer, then the SN is at redshift z = 1.8, and if the latter 
then the redshift is weakly constrained by the pho to-z to 
z = 1. 55 ^ 0(22 ■ After classifying the SN in Sec tion 
will revisit this redshift constraint in Section 14.2 
ing to bear all available information from the 


4.2 


we 


bring- 
'N light 

curve, and also considering the possibility of an unseen 
host galaxy. 
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rest frame time (days from peak) 

-10 0 10 20 30 40 -10 0 10 20 30 40 -10 0 10 20 30 40 



56500 56600 56500 56600 56500 56600 

observer frame time (MJD) 


Figure 7. Infrared light curve of GND13Sto, showing flux as a function of time. Top panels show observations in HST wide-bands, in 
flux units normalized to an AB magnitude zero point of 25. Lower panels show residual fluxes for the corresponding HST medium-bands, 
after subtracting off the predicted flux from the best-fit model. The top axis marks rest-frame time at redshift 2 : = 1.8. relative to the time 
of B band peak brightness. 



F139M - F140W F139M - F140W 


0.4 
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0.2 


0-1 g 
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- 0.2 
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-0.4 


Figure 8. Medium-band pseudo-color circle diagrams for GNDlSSto, as in Figures and In both panels, the observed colors for 
GNDlSSto are indicated by the black square with error bars. (Left) Classification circle diagram: contours enclose 68% and 95% of the 
population for each class, with red for Type la, green for Type H and gold for Type Ib/c. The simulated SNe have redshifts that evenly 
sample the redshift range [1.6,2.4]. (Right) Redshift circle diagram : the line of colored points traces the pseudo-colors of a Type la SN 
from z=1.6 (blue) to z=2.4 (red), with diamond symbols marking each increment of 0.1 in redshift. 
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4.2. GNDlSSto Classification 

In the left panel of Figure we consider the classih- 
cation information provided by the medium-band circle 
diagram (developed in Section]^ and Figure]^. The ob¬ 
served medium-band pseudo-colors of SN GOT)13Sto are 
modestly offset from the origin in this pseudo-color space, 
locating it in the circular region where a Type la SN is 
expected to lie. The uncertainties are large, so although 
Figure]^ is suggestive of a Type la classification, by it¬ 
self this evidence from medium-band pseudo-colors is far 
from definitive. Recall, however, that this information 
from the medium-band imaging is particularly valuable 
because it is insensitive to reddening and it is indicative 
of the presence or absence of specific absorption features 
in the SN’s SED. 

For a more complete classihcation picture, we next 
classify GNDlSSto using a Ba yesian photometric clas¬ 
sification algorithm, as in JR14[ This classiher uses the 
sncosmo software packaga*jto generate simulated SN 
light curves from spec trophotometric t emplates. We rely 
on the SALT2 model (|Guy et al.|2010|) to represent nor¬ 
mal Type la SN, and use 42 spectro photometric tem¬ 
plates from the SNANA template librarjp^ (Kessler et al. 
2009b|) to represent GC SNe (26 TypeU and 16 'I'ype 
Ib/c). Although this template library is at present the 
best available collection of spectrophotometric models for 
CC SNe, it is an incomplete sampling of this hetero¬ 
geneous population. It is possible to incorporate some 
knowledge about this incomplet eness into the photomet - 
ric classihcation framework (e.g. Rodney & Tonry|2009), 
but for this work we make the simplifying assurnption 
that the library provides a nearly complete and evenly 
distributed sampling of CC SN properties. 

For redshifts z > 1.5, the available templates (particu¬ 
larly those for CC SNe) do not extend far enough into the 
rest-frame UV to allow for a comparison to observations 
in ACS-WFC or WFC3-UVIS hlters without extrapo¬ 
lating the templates. Thus, in this analysis we restrict 
our photometric classihcations and light curve htting to 
the available WFC3-IR observations, which correspond 
to rest-frame optical bands. As noted above, we have 
only upper limits for ACS and UVIS fluxes for both of 
these SN, so this stricture does not substantially diminish 
the available information. 

We have made one important change in the classih¬ 
cation algorithm relative to previous work, which is to 
eschew a luminosity prior. To get photometric classih- 
catio ns of the complete CA NDELS/CLASH sample in 
R14 and Graur et al. (2014) we included a soft prior on 
the peak absolute magnitude as a function of redshift 
for each SN sub-class. The prior was based on the ob¬ 
served luminosity functions of these SN types in the local 
universe combined with an assumed cosmological model. 
This was necessary for the calculation of volumetric SN 
rates, as a large fraction of the sample has very limited 
photometric data (in some cases just 1-2 epochs with 2 
or 3 bandpasses). Although in principle one could em¬ 
ploy a weak prior without introducing a signihcant bias, 
in this work we have enough data such that we can get 

sncosmo vl.l, source code at https://github.com/sncosmo 
documentation at https: //sncosmo. readthedocs. org/en/vl .0.x/ 

SNANA: SuperNova ANAlysis package; http://das. sdss2. 
org/ge/sample/sdsssn/SNANA-PUBLIC/ 


robust classihcations without the luminosity/cosmology 
prior. This ensures that these SNe can be used for cos¬ 
mology without any concern of a bias introduced at the 
classihcation stage. 

In the case of GND13Sto, we hrst classify the object 
using both wide- and medium-band photometry but in¬ 
cluding no redshift prior, and hnd that it is strongly clas- 
sihed as a Type la SN. The net classihcation probability 
for all CC SN sub-types is < 1%, and no CC template can 
match the photometric observations with b etter than 
2.8. This agrees with our prior analysis in |R14[ where we 
found that GNDlSSto is robustly classihed as a Type la 
SN based on a Bayesian photometric classihcation using 
the wide-band light curves. 

Our photometric classihcation apparatus can also ac¬ 
commodate a prior that infers the SN clas s based on 
the co lor and morphology of the host galaxy (Foley et al. 
20131. In this case, using either of the plausible host can¬ 
didates, galaxies A and B, we still hnd that GNDlSSto 
is classihed as a Type la SN with > 98% conhdence. 

4.3. GNDlSSto Redshift 

Due to the ambiguity of the host association, it is par¬ 
ticularly valuable to evaluate the redshift constraints on 
GNDlSSto independent of any host. As a hrst step, in 
the right panel of Figure we consider the redshift infor¬ 
mation provided by the medium-band pseudo-color circle 
diagram (see Figure!^. Working now under the assump¬ 
tion that SN GNDlSSto is a Type la SN, the location 
in the upper right corner of this pseudo-color space is 
most consistent with a redshift near z = 1.85 or near 
z = 2.2. Redshifts close to 1.7, 2.0, and 2.4 are disfa¬ 
vored by these data. Due to the low S/N of the medium- 
band observations for this SN, this diagram by itself is 
not dehnitive for redshift determination. However, as we 
will see, the addition of medium band observations to a 
complete wide-band light curve does provide a valuable 
improvement in the redshift precision, and supports an 
association with host galaxy A at z = 1.8. 

In Figurej^we build up a composite redshift constraint 
piece by piece, showing how the posterior probability dis¬ 
tribution over redshift is separately restricted by the host 
galaxy, the wide-band light curve, and the colors from 
a single epoch that includes medium-band observations. 
First, the wide band light curve alone provides a fairly 
weak constraint, z = l.Td^Qg®. Second, from the single 
epoch of medium-band pseudo-colors collected near peak 
brightness we derive a posterior probability distribution 
that is periodic in redshift. With only a single epoch, the 
shape of the light curve is unconstrained (i.e. any change 
in the value of xi could be reproduced by coordinated 
changes in xq and tpk)- To address this, we limit the 
range of tpk to ±3 days around the value inferred from 
the wide-band light curve ht - effectively adding in addi¬ 
tional information by fixing the phase of the light curve 
at the time of the medium-band observations. Formally, 
we derive the redshift and uncertainty as z = 1.90jtg;34, 
but Figure shows that we actually have 4 peaks, at 
z ~ 1.5, 1.9, 2.2, and 2.6. This periodicity arises from 
spectral features moving in and out of the medium bands 
as the SN la template redshift ch anges , illustrated in Fig¬ 
ure and previously observed in S09 By combining the 
wide band light curve and the medium band epoch to- 
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gether, this redshift degeneracy is broken. 

The final piece of redshift information we would like 
to include is a host galaxy redshift prior. Given that 
the host galaxy association is ambiguous, we use a com¬ 
posite redshift constraint that is the sum of the redshift 
probability distributions from host galaxy candidates A 
and B. In effect this allows an equal probability that the 
SN belongs to either galaxy, resulting in a sharp peak 
at z=1.8 (from host candidate A) and broad tails on ei¬ 
ther side (from host candidate B), as shown in Figure]^ 
Including this as a redshift prior, and utilizing all avaiF 
able SN photometry, we get a final redshift constraint 
of ZA-i-B = l-85lo;i2- To test how sensitive this red¬ 
shift is to the host galaxy prior, we have also computed 
the redshift of GNDlSSto using only the photometric 
redshift of host candidate B (l-SSl*!®;®®)- This results in 
a posterior probability distribution that is peaked very 
close to the spectroscopic redshift of host candidate A: 
Zb — l-05_g 

A final possibility to entertain is that the true host 
galaxy of SN GND13 Sto is not detected in the GAN- 
DELS imaging data. Sullivan et al. (2006) found that 
7% of the SNe in the SINLS prograni did not have a 
detected galaxy within Rsn = 5. Setting aside SN 
GNDlSSto the equivalent “hostless” fraction is only 3% 
in the CAND ELS/C LASH programs (Graur et al.||2014 
Rodney et al.||2014D, owing largely to the advantage o: 
very deep HSi' imaging in these survey fields. Never¬ 
theless, to consider the possibility of an unseen host, 
we repeat the redshift test using a flat prior, and find 
-^nohost = 1-951 q; 26- Thus, even if the GNDlSSto host is 
undetected, the SN photometry alone requires the red¬ 
shift to be fairly close to the spectroscopic redshift of 
host candidate A (1.80±0.02). 

We therefore conclude that the most likely scenario 
has SN GNDlSSto as a Type la SN at 2 = 1.80 ± 0.02. 
It may be bound to galaxy A, or associated with some 
possibly unseen component of a group containing that 
galaxy. At redshift z = 1.80 the best-fit SALT2 light 
curve parameters and uncertainties are tpk = 56482.1 ± 
1.2, xi = -0.47 ± 0.68, c = -0.02 ± 0.07, and ms = 
26.14 ± 0.07. Figure!^ shows this maximum likelihood 
Type la model, whicn has = 10-9 for 12 degrees of 
freedom. 

4.4. GNDlSSto Lensing 

The foreground galaxy C at z = 0.5572 is separated 
from the SN by 2/9, making it a possible source of grav¬ 
itational lensing for SN GNDlSSto. Our p rocedure for 
estima ting the lensing magnification follows Jones et al. 

([M^. 

To define the lensing potential, we start with the stel¬ 
lar mass of the galaxy derived from SED fitti ng (Table 
and u se the broken power law relation of [Yang et al. 
(2012) to convert this to the total mass of the dark mat¬ 
ter halo. We then model the stellar com ponent of the 
galaxy with a Sersic profile ( Sersic 1963|) with param¬ 
eters derived using GALFIT (Feng et al. 2002). Next 


we assume a Navarro-Fre nk-White (INT'W) density pro¬ 
file (Navarro et a_L 1997) for t he lensing potential, and 
use GRAYLENS ([Keeton 2001|) to generate 10,000 Monte 
Carlo realizations in which we vary the stellar mass of 
the lensing galaxy using the limits from SED fitting as re¬ 


ported in Table|^ and the redshift of SN GNDlSSto from 
section IL^ Our Monte Carlo simulation also varies the 
mass-concentration relation and t he stellar-to-halo mass 
ratio within their observed ranges (Guo et al. ]MI| |Yang| 


et al.|[20T2 ). From this we infer a median magmhcation 


and 68% confidence limits of 1.02/^g Q^ from the stoneC 
galaxy at the position of SN GNDlSSto, indicating that 
the SN is not substantially biased by gravitational lens¬ 
ing. We include no correction for the insig nificant mag¬ 
nification derived here, but see Section pL^ for discussion 
of lensing effects on a full high-z SN sample. 

5. GND12Col: SN COLFAX 
5.1. GND12Col Host Galaxy 

Our second high-z Type la SN candidate is GND12Col, 
[3 found in the CANDELS GOODS-N field on 27 May 
2012. The nearest galaxy to the SN position, labeled host 
candidate A in Figure has a relatively tight photo-z at 
z = 2.35, with 6 8.3% confidence interval [2.24,2.43]. 
Dahlen et al.| (|^13[) recently demonstrated that even 


well-constrained photometric redshift estimates from any 
single SED fitting algorithm can be systematically biased 
and underestimate the redshift uncertainties, though a 
combination of photo-z estimates from multiple codes 
can reduce these concerns. To address this, we also fit the 
SED of the GND12Col host candidate A using the BPZ 


code, a Bayesian photometric redshift estimator (Benitez 
2000| . We derive a redshift probability distribution tunc- 


tion, P(z), centered at z = 2.17 with 68.3% confidence 
region [2.10,2.23], in good agreement with the photo-z 
derived above. We define the composite photo-z for host 
candidate A as a Gaussian centered at the mean of the 
two, with a 68.3% confidence region encompassing the 
mean of the two photo-z distributions: z = 2.26 ± 0.09. 

A second possible host galaxy - candidate B - is the 
large face-on spiral to the southeast, w'hich has a spectro - 
scopic redshift of z = 0.9466±0.0002 (Gowie et alj2004 ). 


To determine whether galaxy A or B is the more likely 
host, we first evaluate the SN-galaxy separation indepen¬ 
dently of redshift using the R parameter of Equation 
For host candidate A we find R — 1.61, and for host 
candidate B we have R = 4.4, supporting candidate A 
as the more likely host. Including the available redshift 
information strengthens this conclusion. GND12Gol is 
2.8" from the nucleus of host candidate B, which is 23 
kpc at the redshift of that galaxy (assuming a flat AGDM 
cosmology with Ho=70, Di„=0.3). The SN is only 0.44" 
(3.3 kpc at z « 2.2) from the center of host candidate 
A. Based on this evidence, we assign host candidate A 
as the host of SN GND12Col, lending the SN a strong 
redshift prior at z = 2.26 ± 0.09. 

5.2. GND12Gol Glassification and Redshift 


R14 


_ ^classified GND12Col as a Type la SN, and we find 

that this classification holds up with improved photom¬ 
etry and is supported by the addition of medium-band 
observations. The formal Type la classification probabil¬ 
ity is effectively unity, with the combined probability for 
all common GGSN sub-classes limited to < 10“®. Fig¬ 
ure 10 shows the maximum likelihood Type la SN model 


Nicknamed “SN Colfax,” after Schuyler Colfax, 17th Vice 
President of the United States. 
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Wavelength 



Figure 9. The host galaxy of GND12Col. Left: Template F160W image combining all exposures from before MJD=56000 or after 
MJD=56300. Ellipses are drawn at R = 2.5, using parameters measured using SExtractor in the F160W band (see Equation [y. Galaxy 
B is a foreground object. Right: Points show the observed photometry of the host, galaxy A, and the curve shows the best-iftting SEE) 
model. The inset at lower right shows the redshift probability distribution function derived using the LePhare SED fitting code, centered 
at 2 : = 2.26, with a 68% confidence range of [2.17,2.35]. 
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Figure 10. Infrared light curve of GND12Col, showing flux as a function of time. As in Figure [^fluxes are normalized to zero point 25 
and lower panels show residual fluxes for medium-band observations. The top axis marks rest-frame time at redshift 2 : = 2.26. 


at z = 2.26, which has x^=19.8 for 21 degrees of free¬ 
dom. The best-fit SALT2 light curve parameters and 
uncertainties are tpk = 56077.4 ± 5.6, Xi = 0.15 ± 1.06, 
color c = 0.04 ± 0.13, and peak apparent magnitude 
tob = 26.80 ±0.07. 

To visualize how the medium-band observations lead 
to an improvement of the GND12Col classification and 
redshift, we show in Figure [TT] the position of the SN in 
a pair of 2D projections of tHe pseudo-color space. The 
CCSN population is expected to be spread out across 
the entire pseudo-color space, but in both of these slices 
we find GND12Col is conspicuously close to the narrower 
region occupied by Type la SNe. Note also that the SN 
position in both projections of pseudo-color-color space 
is consistent with the middle of the redshift distribution, 
i.e. 0 ^ 2.2. 


Figure [T^shows how the posterior probability distribu¬ 
tion over iTOshift for GND12Gol is separately constrained 
by the host galaxy prior, the wide-band light curve, and 
the colors from a single epoch that includes medium-band 
observations. Each of these components alone provides a 
weak constraint, but all three are consistent and with the 
combination we achieve a precision of =^0 iq, on par 
with the best available photometric redshifts from well- 


samp led galaxy SEDs at similar redshift (Dahlen et al. 


20131. 


5.3. GND12Col Lensing 

As in Section [T^ we evaluate the possibility of gravita¬ 
tional lensing from the foreground galaxy B. This galaxy 
has a projected separation of 278, and a total stellar mass 
of ^10^° Mq. As with GND13Sto, we adopt a Sersic 
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Figure 11. Medium-band pseudo-color diagrams for GND12Col, as in Figures pland]^ Contours enclose 68% and 95% of the population 
for each class, with red for Type la, green for Type II and gold for Type Ib/c. T'he simulated SN have redshifts that evenly sample the 
redshift range [1.9,2.4]. The line of colored points traces the pseudo-color of a Type la SN from z=1.9 (blue) to z=2.4 (red), with diamond 
symbols marking each increment of 0.1 in redshift. Observed colors for GND12Col are indicated by the black square with error bars. 



1.4 1.6 1.8 2.0 2.2 2.4 2.6 

Redshift 

Figure 12. Redshift constraints on GND12Col. As in Figure 
the filled orange curve shows the host galaxy redshift prior (from 
galaxy colfaxA), filled curves in blue and red show the constraints 
from the wide-band light curve and the medium-band colors (with¬ 
out including a redshift prior), and the solid black curve shows the 
composite constraint from all available data. All PDFs are nor¬ 
malized to a peak value of 1.0. 


profile for the stellar component and assume an NFW 
dark matter halo, with mass concentrations and stellar- 


to-halo mass ratios following Guo et al. (2011 ) and Yang 


et al. (2012). From 10,000 Monte Carlo realizations of 
the lensing potential we determine a magnification factor 
ofM = 1.04tH^ Once again we do not include a correc¬ 
tion to the SN magnification for this mild m agnifi cation, 
but examine global lensing effects in Section [7?^ 


6 . EVALUATING THE MEDIUM-BAND 
PSEUDO-SPECTROSCOPY TECHNIQUE 

For the two high-z Type la SNe discussed in this work, 
one may well ask whether the medium-band IR observa¬ 
tions have provided significant value for the SN classifica¬ 
tion and redshift determination. After all, both of these 
objects were classified as Type la SN with high confi¬ 
dence in |R14[ and were assigned redshifts from their host 
galaxies, all without the benefit of the detailed medium- 
band analysis employed here. 

Before scrutinizing the value of the medium band imag¬ 
ing in the following sections, let us first note that relative 
to photometric data the medium-band data are less sus¬ 
ceptible to systematic biases. This is because (a) the 
medium-band pseudo-colors are insensitive to redden¬ 
ing (Figure ® and (b) the redshift information in the 
medium-bana data comes from individual spectral ab¬ 
sorption features (Figure]^, not from the overall shape 
of the SED or the light curve. 


6 .1. Classification 

A key practical benefit of the medium-band approach is 
that it can be used for efficient direction of scarce follow¬ 
up resources toward high-value targets. When a high- 
z Type la SN candidate is initially discovered, a single 
epoch of medium-band imaging can substantially clar¬ 
ify the classification picture while the SN is still close 
to peak brightness - as long as the redshift and time of 
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peak are fairly well constrained. Thus, one can get a 
strong classification assessment for a relatively low cost 
and then evaluate whether the candidate warrants an in¬ 
vestment of additional follow-up time to collect the full 
wide-hand light curves. Importantly, the medium band 
pseudo-colors are insensitive to dust (cf. Figure and 
unaffected by the underlying cosmology. Thus, mSium- 
band observations are similar to spectroscopy in that 
they can be used to select a Type la sample for follow¬ 
up without needing to invoke any priors on luminosity 
or color, which may carry a bias with redshift. The 
medium band classifications are therefore more robust 
than purely broad-band photometric methods, but more 
efficient than spectroscopy. 

In addition to improving the SN redshifts, the medium- 
band data also add value by strengthening our photo¬ 
metric SN classifications. To see how large this effect 
is, we have done a piecemeal classification test: running 
subsets of each SN light curve through our Bayesian pho¬ 
tometric classification algorithm both with and without 
the medium-band imaging data. Figure [T3| shows the re¬ 
sults. The value of the medium-band data in strengthen¬ 
ing the classification is apparent in the first two columns 
of this figure. For these two tests we used just a single 
wide-band light curve (FI25W or FI60W) plus a sin¬ 
gle epoch of medium-band imaging (including the nec¬ 
essary observations in each complementary wide band). 
The addition of medium-band imaging alters the prob¬ 
abilities from a misleading or ambiguous P(Ia|D)<0.5 
into a much more robust P(Ia|D)>0.75. In these cases 
the single-band light curve primarily constrains the time 
of maximum light, while the single-epoch medium-band 
pseudo-colors provides the majority of the classification 
information. In the third column of Figure |I3[ where the 
full multi-band IR light curve is available, the medium- 
band pseudo-colors do not alter the formal classification 
probability. This shows that the Type la classification 
can be securely established without any medium-band 
data, but to do so requires shape and color information 
afforded by multiple wide-band light curves extending 
over several months. 

6 .2. Redshift Determination 

In Sections 14.31 and 15.21 we saw that the medium-band 
imaging improves the aecuracy of the final redshift de¬ 
termination, by providing redshift constraints that are 
independent of the presumed host galaxy. In the case of 
SN GNDISSto, where the SN is widely separated from all 
nearby galaxies, the additional redshift constraint from 
medium-band pseudo-colors is especially valuable as a 
check against mis-identification of the host galaxy. For 
SN GNDI2Col “ where we have no spectroscopic data on 
the host galaxy - the redshift information derived from 
medium-band imaging is very valuable to guard against 
a catastrophic error in the host galaxy photo-z, or misas- 
sociation with a galaxy that is actually in the foreground 
or background. 

Figures [6| and also demonstrate the value of medium 
band imaging to improve the precision of redshift con¬ 
straints. For both GNDISSto and GNDI2Col the 
medium-band pseudo-colors deliver posterior redshift 
distributions with multiple peaks. As the broad single 
peak from the wide-band light curve breaks the degener¬ 
acy by picking out a single peak, the combination of wide- 


FI60W FI25W All Wide 
Light Curve Light Curve Band IR 



Figure 13. Examining the role of medium-band data in modify¬ 
ing the Type la classification probability for SN GNDISSto (blue 
diamonds) and GND12Col (red squares). Along each column we 
plot the Type la classification probability derived from a subset 
of the light curve data. For column 1 we use just the F160W 
band light curve (with and without a single epoch of medium- 
band pseudo-colors). Golumn 2 uses just the F125W band light 
curve, and column 3 uses all available wide-band IR data. Filled 
symbols with solid lines mark the results when including medium- 
band imaging data, and open symbols with dashed lines show the 
probabilities inferred without any medium-band data. In all tests 
we include a redshift prior derived from the nearest host galaxy 
candidate, but we add no prior on the time of peak brightness. 


and medium-band imaging can reduce the photometric 
redshift uncertainty by almost a factor of 2. Redshifts 
measured from the host galaxy can improve the precision 
still further, but if either of these SNe had no apparent 
host galaxy, then the medium band imaging would have 
been critical to make the SN useful for cosmology. 


6.3. Future Surveys 


For SN GNDISSto we collected exposures totaling ~5 
ksec (2 HST orbits) for each of 2 medium band filters 
(FI39M and FI53M). For GNDI2Gol we used roughly 8 
ksec (3 HST orbits) on each of 3 medium-band filters. In 
both cases we achieved a final S/N ratio between 5 and 
10 in each individual filter, and we have seen that this is 


insufficient t o p rov ide a definitive classification by itself 
(e.g. FiguresjMand II). To evaluate medium-band imag¬ 
ing for future surveys, it is instructive to consider how 
much exposure time would have been required to make 
the medium-band imaging pseudo-colors stand alone and 
deliver a robust classification without reference to the 
wide-band light curve at all. 

To that end, we performed a Monte Garlo classification 
test in which we simulated medium-band observations for 
2500 synthetic SNe spread evenly across the primary SN 
sub-classes, with redshifts between 1.6 and 2.4. Each 
simulated SN was “observed” through 3 medium bands 
and their complementary wide-band filters, with all ob¬ 
servations taken concurrently in a single epoch within 
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Figure 14. Results of a Monte Carlo simulation test examining the use of medium-band pseudo-colors alone for SN classification. 
Histograms show the distribution of Type la classification probabilities derived by classifying 2500 synthetic SN ’’test particles” using only 
a medium-band pseudo-colors. All SNe are simulated at redshifts in the range [1.6,2.4]. A single epoch of simulated observations are 
made in three medium band filters (F127M, F139M, F153M) and their corresponding wide-band counterparts (F125W, F140W, F160W), 
taken at ±15 observer-frame days from peak brightness. The blue histograms show classification results for simulated CC SNe and the red 
histograms count probabilities for synthetic Type la SNe. The test particles were simulated in three rounds, with the flux uncertainty in 
every band fixed to a S/N ratio of 10, 20, or 30 — shown in three panels from left to right. 


±15 observer-frame days of peak brightness. We then 
derive a Bayesian classification probability P(Ia|D), tak¬ 
ing care to avoid self-identification, where a synthetic SN 
yields an artificially good match when compared with the 
template that was used to create it. For CC SN test parti¬ 
cles, this means that we remove from the template library 
the single template used to generate the CC SN test par¬ 
ticle. For Type la SN test particles, we exclude from our 
comparison the range of SALT2 light curve shape and 
color parameters that are similar to the values generated 
for the test object. Here, “similar” is arbitrarily defined 
as a circle in (xl, c) phase space centered on the test 
particle’s (xl, c) values, with radius 0.05. 

Figure [l^ shows histograms of the inferred type la 
classification probabilities, derived from a series of three 
Monte Carlo simulations. In each simulation we fixed the 
S/N to be 10, 20 or 30 for every individual filter. The 
progression of three panels in Figure 14 demonstrates 
the result that is to be naturally expected: as S/N is 
increased, the prospect of classification from medium- 
band data alone becomes more tenable. With S/N=10 
in each band, the two simulated populations are indis¬ 
tinguishable. At S/N=20 the majority of synthetic Type 
la SNe are recovered, with P(Ia|D)>0.75 - but there is 
still substantial contamination from CC SNe being incor¬ 
rectly classified. A large Type la sample defined using 
medium-band observations at S/N=20 would have a pu¬ 
rity of ~80%. When all observations attain S/N=30, the 
populations are much more distinct, and it would be pos¬ 
sible to define a Type la SN sample with >95% purity 
using such high-quality medium-band data. 

Realistically, these S/N targets are not achievable with 
HST for targets like SN GND13Sto and GND12Col. To 
reach S/N=20 for a Type la SN at z=2.0 with an F160W 
AB mag of 26 would require a total of roughly 100 ksec 
(40 HST orbits) spread across three medium bands and 
their wide-band counterparts. For S/N=30, the expo¬ 


sure time requirement rises to ~215 ksec (nearly 90 HST 
orbits). At those costs, grism spectroscopy may become 
an appealing alternative, especially since it delivers more 
ancillary science benefits in any given field. Given the 
limitations of finite observing time, the best strategy for 
HST medium-band imaging may be the one used in this 
work: relatively short observations that are not definitive 
in themselves, but provide support for a holistic approach 
to photometric classification and redshift determination. 

The small set of medium-band IR filters on HST cur¬ 
rently limits the application of the medium-band pseudo¬ 
spectroscopy method to SNe at z ^ 2. The JWST 
Near-infrared Camera (NIRCAM) will have 12 medium- 
band filters from 1.5 to 5 /rm, substantially expanding 
the reach of this technique and the potential for testing 
models of time variable dark energy or Type la SN evo¬ 
lution. In principle, JWST medium band filters could 
be used to characterize Type la SN candidates as far as 
z 6. Unfortunately, the diminishing cosmic star forma¬ 
tion rate density at very high redshifts will make these 
objects few and far between, so the initial discovery will 
be more efficiently done with a wide-held observatory, 
such as Euclid or WFIRST-AFTA. With this combina¬ 
tion of Euclid/WFIRST and JWST, it will be possible 
to discover and characterize a sample of primordial Type 
la SNe- explosions of the very hrst generation of binary 
white dwarfs - which in turn would enable the most strin¬ 
gent possible test for the secular evolution of Type la SN 
properties. 

7. THE HIGH-Z FRONTIER 
FOR TYPE lA SN COSMOLOGY 

The two objects presented here are among the most 
distant Type la SNe known, with GND12Col now set¬ 
ting the bar with the highest redshift. As such, these 
SNe are of great interest for testing cosmological models 
and evaluating possible systematic biases that may be 
suppressed or obscured at lower redshifts. In this sec- 
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tion we b egin with an extension of the Hubble diagram 
in Section |7.H and then examine some of the particular 
considerations that must be addressed when incorporat¬ 
ing such high-z SNe into a cosmological sample. 

7.1. Extending the Hubble Diagram 

In Figure[^we plot a Hubble residuals diagram, show¬ 
ing the luminosity distance relative to ACDM versus red- 
shift. The baseline cosmology is defined by fittin g to a 
collection of 943 Type la SNe from the literature ([Riess 
et al. 1119991 |Jha et al.||2006| IHicken et al.||2009| IKessler 
-- - 'lli et al 


et al. 2009a 


Contreras et al.||2010[ |Folatelli et al.||2010[ 


Conley et al. 20111 Stritzinger et al. 20111 Hicken et ^ 


2012 Caneshalingam et al. 12013 IRest et al. 2014|). We 

fit these light curves using the SAL'i'2 model, and define 


a corr e cted a pparent B-band peak magnitude following 
Tripp (19981: = tob + axi — j3c. Here the (un- 

correctedj peak apparent magnitude me, the light curve 
width xi, and the color c are derived from the light curve 
fit. For this exercise we fixed t he SALT2 model parame¬ 
ters to the values derived in the Rest et al. (2014) analysis 
of the PSl SN sample: a = 0.147 and p = 3.13. We then 
subtract mj^cDM^ the apparent magnitude predicted by 
a flat ACDM cosmology at each redshift, with the Type 
la SN peak absolute magnitude at a fiducial value of 
Mb = —19.34, and using best-fit values of Hm=0.284 
and iJo=70.88. The SN magnitude errors shown in Fig¬ 
ure also include an intrinsic scatter of a\„t = 0.12 
magnitudes ( Scolnic et al.|2014 Betoule et al.|2014 |, and 
for simplicity this term is not separately defined tor each 
survey. These simplifications allow for an expedient com¬ 
parison of the new high-z Type la SNe against existing 
SN and a baseline cosmological model. A more complete 
examination of cosmological models will be presented in 
a subsequent paper (Riess et ah, in preparation). 

7.2. Gravitational Lensing 

As Type la SN observations are pushed to higher red- 
shifts, observed SN fluxes are increasingly affected by 
gravitational lensing. Photons in the Universe are con¬ 
served, so when sampling a sufficiently large number of 
random sight lines, the mean SN flux is unperturbed. 
Thus, lensing adds scatter to Type la SN distances, but 
if treated appropriat ely it does not bias cosmological in¬ 
ferences ( ]Hoiz||1998 |. 

Lensing can be handled in one of two ways: either (a) 
rely on a large sample to average out the bias, while ac¬ 
counting for the increased scatter, or (b) estimate the 
(de)magnification along every line of s i ght a nd include 
appropriate corrections. Holz & Linder (2005) examined 
the former option and found that lensing introduces an 
additional redshift-dependent dispersion in the observed 
Type la SN fluxes of ciensing = 0.088z (or 0.093z if treat¬ 
ing intrinsic Type la SN dispersion as Gaussian in mag¬ 
nitude). 

A formalism for i ndividual lensing correctio ns (option 
b) was presented in Gunnarsson et al. (2006), in which 
galaxy detections from deep imaging surveys are used 
to mo del dark matter halos that intersect each SN sight 
line. Jonsson et al. (|2006|) applied this method to 33 
SNe from the HS'l ' GUUDS survey ( Strolger et al.|[2004 


ditionally, [Jonsson et al. found that the mean magni¬ 
fication factor in the HST sample was consistent with 
unity, implying that there is no evidence for a bias to¬ 
ward including highly magnified objects in the high-z SN 
cosmology sample. 

In this work we have found no evidence for significant 
lensing from the n eare st de tect ed foreground galaxy to 
each SN (Sections 4.4 and 5.3). However, the nearest 
foreground galaxy is not the only possible source for grav¬ 
itational lensing perturbations. For a more precise lens¬ 
ing correct ion, one should include all foreground source s 
out to ^l'( Gunnarsson et al.||2006 Jonsson et al.||2006 ). 

As we are not pursuing a measurement of cosmological 
parameters in this work, we adopt option (a) from a bove, 
applying t he 0.0 93z mag additional scatter from Holzj 
& Linder (2005) as an additional noise term for all SN 


Riess et al. 2004b) and demonstrated that such iensmg 
corrections can indeed reduce the lensing scatter. Ad- 


shown in fi igure|l5| For our two high-z Type la SNe this 
corresponds to 0.17 mag at the redshift of GND13Sto 
and 0.2 mag for GND12Gol. It is noteworthy that in 
this high redshift regime we have transitioned to a point 
where the scatter from gravitational lensing along the 
line of sight is larger than the intrinsic scatter of the 
Type la SN population (~ 0.12 mag)-though for these 
individual objects both of those error terms are dwarfed 
by the observational uncertainties. 

7.3. Biases and Evolution 

At the high-z edge of any survey, a selection bias 
may be introduced as intrinsically faint objects are more 
likely to slip below the detection threshold and be missed 
( Malmquist|1920 ). The blue shaded regions in the upper 
right corner of 1 igure 15 show the redshift and luminos¬ 
ity regime where such abias could begin to remove Type 
la SNe from our sample. These regions are defined by 
simulating a fiducial Type la SN {xi = 0, c = 0) at each 
redshift and measuring the visibility time, tvis- the num¬ 
ber of observer-frame days where the SN is brighter than 
the 50% detection efficiency threshold of the GANDELS 
SN search (26.5 AB mag in F125W and F160w, seejR14 ). 
In the light blue region, our fiducial Type la SN model 
has a visibility time shorter than the typical cadence of 
the GANDELS survey (50 days). In the dark blue region, 
the visibility time drops below zero, meaning that normal 
Type la SNe are effectively undetectable. This does not 
paint a complete picture of the Malmquist bias, as these 
regions depend on the light curve shape and color, but 
this illustration shows that the absence of positive Hub¬ 
ble residuals at z > 1.5 should be unsurprising, given the 
measured detection limits of this survey. 

In a subsequent paper, we will present distances for the 
complete sample of high-z Type la SNe from GANDELS 
and GLASH, comprising ~15 Type la SNe at 1 < z < 1.7 
(Riess et al. in preparation). Further high-z Type la SNe 
will be delivered by ongoing HST surveys such as the 
FrontierSN program (PLRodney) and the See Change 
survey (PLPerlmutter). These extensions of the Type la 
SN sample at z ^ 1.5 will improve the sensitivity of Type 
la SN cosmology, especially in testing for a time-variable 
dark energy equati on of state (the ratio of pressure to 
density, w = P/p). jSalzano et al.j (|2013|) quantified how 
much these and future surveys can be expected to refine 
our constraints on w. Adopting a common parameteriza- 
tion that is linear with the scale factor, w = Wn-|-Wn.(l—a) 

( Chevallier fc Polarski]200T Linder]2003 ), Salzano et al. 
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Figure 15. Hubble residuals diagram, plotting as a function of redshift the apparent magnitude (corrected for light curve shape and color) 
relative to the predict ed magnitude from a flat ACDM cosmolog y. Colo red circles show a co mpilat ion of ^900 SNe from th e low-redshift 
CfA a nd CSP surveys (|Hicken et al.j2009||Stritzinger et aL|2011|l, SDSS ||Kessler et al.j2009a||, PSl jRest et al.|2 014ll' SNLS jConley et al. 
2^^Rj_and at z > 1 the orange circles are trom the UJUUJJb and SCP surveys using me ACS camera on tilij r l|Hiess et al.|2UU7Hbuzuki 
et al.||2012[|. Red diam onds at z=1.55 and 1.91 plot two previously published Type la SNe from the CANDEls survey dRodney et aE 
:^Ui:^|]jones et al.||2013|l, and dark red diamonds at z = 1.8 and 2.26 show the two new Type la SNe from this work. All vertical error 
bars include a unirorm intrinsic dispersion of crjnt = 0.12 mag and a redshift-dependent gravitational lensing dispersion of cr^L = 0.093^:. 
For SN GND13Sto and GND12Col the x-axis error bar is shown as a curved line to indicate the covariance of distance with redshift for 
this object (dashed for GND13Sto to indicate that this redshift range only applies if the association with host galaxy A is presumed to be 
spurious). The shaded region in the upper right demarcates the space where the CANDELS survey starts to lose sensitivity for discovery 
of Type la SNe (see text for details). 


found that a sample of ~ 30 Type la SNe at z > 1 should 
reduce the uncertainty on Wa by ~ 20%. 

It is naturally expected that the demographics of the 
observed Type la SN population should change with red¬ 
shift, as brighter Type la SNe with slower light curves are 
more prevalent in star -forming galaxies, w hich dominate 
the universe at z > 1 Howell et al. (2007). Such a shift 
in the composition of the 'i'ype la SbJ population should 
not by itself introduce any cosmological bias. However, 
it is also possible that differences in the progenitor sys¬ 
tems for high z Type la SNe could lead to luminosity 
biases - deviations from the expected standard candle 
relations between shape, color and luminosity. Indeed, 
one of the original motivations for this high redshift SN 
observing program was to detect SNe in the early uni¬ 
verse where such e ffects could be disen tangled from cos¬ 


mological changes ( jRiess fc Livio|2006 |. 

U sing t he full CAiNUULS and C ^ASH Type la SN sam¬ 
ple, R14| found evidence that the fraction of Type la SNe 
exploding within 500 Myr after formation is < 50%. This 
means that a substantial fraction of all the Type la SNe 
that we might observe at z ^ 2 could have been formed 
when the universe was < 2 Gyr old, perhaps in very low 
metallicity environments. The mean mass of a Type la 
SN progenitor star must also increase with redshift, as 
slow-evolving low-mass progenitors have less time avail¬ 
able to leave the main sequence and start accreting mass 


as a white dwarf. It is conceivable that such changes in 
the progenitor mass or metallicity could lead to a sys¬ 
tematic shift in the peak magnitudes that is not cor¬ 
related with the light curve width or color (iDomlngue z 

llonni 1 .^1 llonnol r 


2003 Riess &: Livio' 200fc)p . 


et al. 2001 Timmes et al. 

F inding evidence for this eflect is more challenging than 
improving constrain ts on variable dark energy models. 


Salzano et al. (2013) found that testing plausible models 
for non-cosmological evolution of the Type la SN popula¬ 
tion requires ~50 SNe at 1.5 < z < 3.5, which could first 
be done with the James Webb Space Telescope (JWST), 
scheduled for launch in 2018. 


8 . SUMMARY 

In this work we have presented two primary results, 
where the first is methodological and the second obser¬ 
vational: 


1. a description of the use of HST medium-band IR 
filters for analyzing high-z SNe; and 

2 . the discovery, classification, and redshift determi¬ 
nation for two of the most distant Type la SNe yet 
known. 

We have demonstrated that the medium-band tech¬ 
nique can provide an alternative to costly spectroscopic 
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observations of high-z SNe. This method capitalizes on 
the semi-periodic nature of the broad absorption fea¬ 
tures in Type la SN SEDs and provides a photometric 
tool for examining spectroscopic features. For classifying 
and determining redshifts of high-z SNe, this technique 
is more efficient than spectroscopy but more precise than 
photometry. By relying on difference imaging, we avoid 
contamination from host galaxy light that can plague 
SN spectroscopy. By using medium-minus-wide pseudo¬ 
colors, this approach does not sample the broad-band 
SED colors and is therefore insensitive to dust extinc¬ 
tion. 

Eor the second component of this work, we have mar¬ 
shaled all available wide- and medium-band photome¬ 
try and host galaxy information to determine that SN 
GNDlSSto and GND12Col are both Type la SNe at 
z ~ 2. These objects push the frontier for Type la SN 
cosmology into an era when the Universe was only 3 Gyr 
old. Although the uncertainties are large, we find that 
both are fully consistent with a standard flat AGDM cos¬ 
mology, showing no obvious deviations that might arise 
from extreme evolution of the Type la SN progenitor 
population. We leave a more robust examination of cos¬ 
mological models to future work, which will incorporate 
the full GANDELS and GLASH SN samples (Riess et al. 
in preparation). 

Whereas the SN samples at z < 1 are moving quickly 
toward the “big data” era, the high-z regime is still very 
much limited by the sample size. With these high-z ob¬ 
jects at the redshift frontier, each requires careful atten¬ 
tion to accommodate the unique nuances of classification 
and redshift estimation. Infrared imaging surveys from 
HST, JWST and successor observatories should continue 
to expand the sample of such high-z SNe in coming years. 
A decade from now the high-z regime may well look like 
the z < 1 samples of today, with hundreds or thousands 
of well-sampled Type la SNe from observatories such as 
NASA’s planned WFIRST-AETA. That will open up the 
exciting possibility of exploring a wider range of cosmo¬ 
logical models and evolution scenarios with these early- 


universe SNe. 
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APPENDIX 

GSDllBus 


Detected in the GANDELS GOODS-S field, GSDllBuj^ appeared in a very low surface brightness host galaxy 
(26.8 AB mag arcsec“^ in the F160W band). The host also nas a red SED, with a V-H color of approximately 1 mag. 
Identification of the 4000 A break across the riz bands leads to a central peak in the photo-z at z = 1.15 as shown in 
Figure [THl However, the large photometric uncertainties allow for other possible solutions at z ~ 0.3, 0.8 and 2.1. No 
spectroscopic constraint on the host redshift is available. 

The SN was initially classified as a possible SN la at z > 1 based on its very red optical-IR colors. This prompted 
follow-up imaging with HST, using both medium and wide bands. However, with the full l ight c urve this object was 
found to be inconsistent with being a normal SN la anywhere in the allowed redshift range. |R14| classified GSDllBus 
as a Type Ib/c, with a redshift of z = 1.7 ± 0.1 based on the combination of the host photo-z and the SN light curve pH 
Using our improved photometry (Tablel^, we revisit this classihcation, and find again that GSDllBus is incompatible 
with SN la models. The classification prwabilities are P(Ia|D)=0.03, P(Ib/c|D)=0.92, and P(H|D)=0.05. The best-fit 
light curve m odel shown in Figure [l7| is based on the Type Ic SN 2006fo from the SDSS SN survey (SDSS ID 13195, 


Zheng et al.|2008). At z = 1.15 this model provides a good ht to the optical-PlR light curve, with XZ = 0.85. 


I'igure |18[shows the redshift constraints on GSDllBus from the wide-band light curve, medium-band colors, the 
host photo-z, and the combined constraint from all sources. Without the host prior, the light curve constraints would 


Nicknamed “SN Bush” after the 41“* and 43’^'^ Presidents of 
the Uni ted St ates. 

In |R14| the photo-z for the GSDllBus host galaxy was re¬ 


ported as z = 1.76 di 0.53, based on an earlier version of the SED. 
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SN GSDllBus 

03:32:42.776, -27:48:07.1 

F160W 

Template 
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Figure 16. The host galaxy of GSDllBus. Left: The F160W template image, constructed from all exposures before MJD=55720 and 
after 55600. The ellipse encloses ~90% of the host galaxy flux in the F160W band. Right: Filled points show the observed SED of the 
GSDllBus host galaxy and the gray curve is the best-fitting SED model. The inset shows the photometric redshift probability distribution, 
peaked at z=1.15. 


Table 5 

GSDllBus Photometry 


Obs. Date 
(MJD) 

Filter 

Exp. Time 
(sec) 

Flux 

(counts/sec) 

Flux Err 
(counts/sec) 

AB Mag^ 

Mag Err 

AB Zero Point 

AZP*’ 

(Vega-AB) 

55720.1 

F606W 

5211 

0.0623 

0.0497 

29.506 

0.865 

26.493 

-0.086 

55723.0 

F606W 

5211 

-0.2223 

0.0546 

>28.5 


26.493 

-0.086 

55795.2 

F606W 

1007 

0.0968 

0.1053 

>27.7 


26.493 

-0.086 

55720.2 

F814W 

5808 

-0.0921 

0.0633 

>27.8 


25.947 

-0.424 

55770.9 

F814W 

1886 

0.3049 

0.0864 

27.237 

0.308 

25.947 

-0.424 

55795.1 

F814W 

932 

0.4176 

0.1057 

26.895 

0.275 

25.947 

-0.424 

55847.7 

F814W 

480 

-0.4589 

0.1879 

>26.6 


25.947 

-0.424 

55720.0 

F850LP 

2046 

-0.0480 

0.0566 

>26.8 


24.857 

-0.519 

55722.9 

F850LP 

2046 

-0.2579 

0.0694 

>26.6 


24.857 

-0.519 

55803.2 

F098M 

2512 

0.4756 

0.0766 

26.474 

0.175 

25.667 

-0.562 

55803.9 

F105W 

1256 

1.2029 

0.1465 

26.068 

0.132 

26.269 

-0.645 

55892.5 

F105W 

8395 

0.2638 

0.0774 

27.716 

0.319 

26.269 

-0.645 

55718.1 

F125W 

1006 

0.1360 

0.1691 

>27.0 


26.230 

-0.901 

55778.8 

F125W 

1006 

0.9312 

0.1593 

26.308 

0.186 

26.230 

-0.901 

55795.1 

F125W 

1006 

1.0172 

0.1556 

26.212 

0.166 

26.230 

-0.901 

55803.1 

F125W 

1256 

1.2903 

0.1558 

25.953 

0.131 

26.230 

-0.901 

55819.6 

F125W 

1006 

0.9515 

0.1565 

26.284 

0.179 

26.230 

-0.901 

55854.2 

F125W 

1006 

0.5491 

0.1548 

26.881 

0.306 

26.230 

-0.901 

55869.9 

F125W 

1006 

0.3475 

0.1442 

27.378 

0.451 

26.230 

-0.901 

55923.2 

F125W 

1006 

0.4099 

0.1440 

27.198 

0.382 

26.230 

-0.901 

55974.6 

F125W 

1006 

0.2356 

0.1477 

27.800 

0.681 

26.230 

-0.901 

55803.3 

F127M 

2512 

0.2966 

0.0658 

25.961 

0.241 

24.641 

-0.961 

55803.3 

F139M 

2512 

0.3664 

0.0725 

25.569 

0.215 

24.479 

-1.079 

55803.9 

F140W 

1256 

1.6063 

0.2136 

25.937 

0.144 

26.452 

-1.076 

55847.8 

F140W 

812 

0.9162 

0.2536 

26.547 

0.300 

26.452 

-1.076 

55803.8 

F153M 

2512 

0.3983 

0.0776 

25.463 

0.211 

24.463 

-1.254 

55718.1 

F160W 

1056 

0.0097 

0.1468 

>26.8 


25.946 

-1.251 

55778.8 

F160W 

1006 

0.7075 

0.1900 

26.322 

0.291 

25.946 

-1.251 

55795.1 

F160W 

1006 

1.0666 

0.1808 

25.876 

0.184 

25.946 

-1.251 

55803.1 

F160W 

1256 

1.0519 

0.1437 

25.891 

0.148 

25.946 

-1.251 

55819.6 

F160W 

503 

0.8344 

0.2442 

26.143 

0.318 

25.946 

-1.251 

55854.2 

F160W 

1156 

0.9717 

0.1809 

25.977 

0.202 

25.946 

-1.251 

55869.9 

F160W 

1056 

0.5799 

0.1598 

26.538 

0.299 

25.946 

-1.251 

55923.2 

F160W 

1006 

0.0633 

0.3185 

>26.0 


25.946 

-1.251 

55974.6 

F160W 

1206 

0.1534 

0.1777 

>26.6 


25.946 

-1.251 


^ For flux values of less than Icr significance we report the magnitude as a 3cr upper limit 
^ Zero point difference: the magnitude shift for conversion from AB to Vega magnitude units. 
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Figure 17. Infrared light curve of GSDllBus, showing flux as a function of time. As in Figures [ t] and [To] fluxes are normalized to zero 
point 25 and lower panels show residual fluxes for medium-band observations. The top axis marks rest-frame time at redshift z = 1.15. 
Grey curves show the best-fit model, based on the Type Ic SN 2006fo. 



Figure 18. Redshift constraints on GSDllBus. As in Figuresand [T^ the independent constraint from each of 3 different sources and 
the combined constraint are shown as probability distribution functions over redshift, all normalized to a peak value of 1.0. 


favor a redshift at z = 1.6 or 0.8. These redshifts also require a CCSN classihcation, so although the host and light 
curve constraints are not in agreement, there is no plausible solution in which GSDllBus can be classified as a Type 
la. The Type la classification is strongly disfavored principally because the slow decline of the light curve favors a 
high redshift, which is incompatible with the observed fluxes in the ACS F814W and WFC3-1R F105W bands. This is 
reinforced by the medium-band colors, shown in Figure [T^ which locate GSDllBus in a region of pseudo-color space 
that is well removed from the Type la SN locus. 

Although this evidence is sufficient to rule out a Type la classification, one may be concerned that the GSDllBus 
medium band pseudo-colors place it outside of the 95% contours for all SN classes, as shown in the left panel of 
Figure [T^ It is possible that GSDllBus is simply a statistical outlier, but we must caution that this discrepancy may 
be indicative of a fundamental problem with the library of CC SN templates employed in this work. Alt hough the 


et al. 

1999 Di Garlo et al. 

2002 

Baron et al. 

2004 

). Fach template was generated by warping those tew spectral 

time series to match the we 

1 -sampled broad-band p 

notometry of 26 low-z SN. In this way, the broad-band colors of 


our simulated SNe can be expected to accurately represent the diversity of broad-band colors from this population. 
However, the diversity of emission line strengths may not be well represented in this library, as the warping procedure 
is not sensitive to such features. 
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Figure 19. Medium-band color-color diagrams for SN GSDllBus, as in Figure!^ Contours enclose 68% and 95% of the population for 
each class, with red for Type la, green for Type II and gold for Type Ib/c. The sirnmated SN have redshifts that evenly sample the redshift 
range [0.7,1.5], encompassing the two primary peaks of the host galaxy photo-z. Observed colors for GSDllBus are indicated by the black 
diamond with error bars. 
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Figure 20. Demonstrating the lack of spectral diversity in the Type II template library employed in this work. All Type II spectrophoto- 
metric templates are shown as blue lines, redshifted to z = 1.35 and observed at +6 days after peak brightness (a redshift and age consistent 
with the SN GSDllBus data). Each template spectrum is normalized to have the same integrated flux through the F160W bandpass, 
which is shown as a solid red line. At this redshift the Ha line occupies the entire F153M bandpass, which is shown as the orange shaded 
region, but our library of models exhibits relatively little variation of the Ha equivalent width. 
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Although we have classified GSDllBus as a Type Ib/c object, the light curve may also be consistent with the Type 
II-L or Iln sub-classes, which are under-represented in our template library. For those cases, the strong Ha feature 
seen in most Type II SNe may be particularly relevant to SN GSDllBus, as that feature intersects the F153M band 
at redshifts near z = 1.35, close to the peak of the GSDllBus host photo-z distribution. Figure shows that our 
template library has a marked lack of diversity in Ha flux relative to the integrated flux through the encompassing 
broad band. Direct observations have shown a much wider range of diversity in the strength and shape of the Ha 


feature (|Patat et a l.|19 94 Gutierrez et al.|2014|), as would be expected for this heterogeneous class. This suggests that 
pnc 


alp 

[JSU 


the medium-band photometry for SIN GSDllBus could lead to a more definitive classification and redshift if we had 
a better set of GG SN templates that more accurately represents the true diversity of these SN types. 
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